Entropic 13 Moment Equations from Boltzmann's Kinetic Equation 
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We use guiding principles from nonequilibrium thermodynamics to develop an admissible set of 
13 moment equations for rarefied gas flows. The main benefits of our thermodynamic approach are 
an explicit entropy expression fulfilling an H theorem and a sound Hamiltonian formulation of the 
reversible free flight transport. To calculate the entropy and to find explicit closure approximations, 
we propose a simple set of approximate 13 parameter solutions to Boltzmann's kinetic equation. 
We discuss how standard hydrodynamics is recovered as a limiting case. 



Introduction. — Rarefied gas flows play an important 
role in a variety of applications ranging from microflu- 
idics to aerodynamics of satellites and space stations in 
the outer atmosphere. For these applications, the Navier- 
Stokes-Fourier equations of hydrodynamics are insuffi- 
cient and the use of Boltzmann's kinetic equation is 
computationally too time consuming. Moment equations 
hence offer an attractive intermediate level of description 
for rarefied gas flow. 

The idea of moment equations goes back to Grad's pi- 
oneering work [2] on solving Boltzmann's kinetic equa- 
tion. For a description of both momentum and heat flux, 
one usually employs 13 moments. In recent years, a prac- 
tical, stable, and accurate set of regularized 13 moment 
equations has been developed by combining the ideas of 
the Chapman-Enskog and Grad methods (see, for ex- 
ample, the textbook Q, the letter Q, and references 
therein). However, several fundamental issues are still 
open, most prominently, the existence of an entropy be- 
yond the linear case and the hyperbolicity of the descrip- 
tion of free flight transport [5]. Closure based on max- 
imum entropy distributions has been suggested to solve 
these remaining problems @, but this idea seems to be 
limited to the 10 moment case, that is, to the absence of 
heat flow @,S). 

The purpose of this work is to develop a fully nonlinear 
set of 13 moment equations with the help of nonequilib- 
rium thermodynamics f§-filj|. We here use two guiding 
principles: (i) the existence of an entropy conserved by 
reversible processes, and (ii) the possibility of a Hamil- 
tonian formulation of reversible dynamics. For our final 
equations, we obtain an entropy satisfying an H theorem 
and a proper mechanical description of the reversible evo- 
lution associated with free flights. 

Variables. — We start from Boltzmann's kinetic equa- 
tion for the single-particle distribution function / = 
f(r,p), that is, the probability density for finding a par- 
ticle with momentum p at the position r. The normal- 
ization is such that n(r) — J f(r,p)d 3 p is the particle 
number density. Angular brackets (. . .) denote position 
dependent momentum space averages performed with the 
normalized probability density f /n. 



We here choose the mass density p = nm (where m is 
the particle mass), the momentum density M — nip), 
and the symmetric second moment tensor tt = (pp) — 
(p) (p) as independent fields. From the third moment 
tensor of rank three, Q = ((p — (p))(p — (p))(p — (p))), 
we construct the vector q = Q : tt^ 1 to arrive at our 
complete set of 13 fields (p, M, tt, q) for which we seek 
evolution equations. Note that the third-rank tensor Q 
is double contracted with tt' 1 rather than the unit ten- 
sor 1. Whereas the latter option would be motivated 
by the form of the heat flux vector, our choice of q 
makes the Hamiltonian formulation of reversible dynam- 
ics more transparent. The use of general vectors and 
tensors (rather than Cartesian ones), with the possibil- 
ity of distinguishing between contravariant and covariant 
objects, is known to be essential for a deep analysis of de- 
formation and flow [L2\ . Just like M, the variable q is 
introduced as a covariant general vector, and tt is a co- 
variant tensor. 

Energy and entropy. — For our two guiding principles 
from nonequilibrium thermodynamics, energy and en- 
tropy play important roles. In terms of our independent 
system variables, the total kinetic energy of all nonintcr- 
acting particles in the volume V occupied by our rarefied 
gas is given by 



E = 



M 2 
~2p~ 



2m 2 



tr tt d A r. 



(1) 



For the nonequilibrium entropy of the ideal gas on the 13 
moment level, we assume the general form 




2^\ 3 m 2 , 
ttt — det TT 



h 2 J 



+ 2 + 



pd 3 7 



where &)b and h are Boltzmann's and Planck's constants, 
respectively, and <p = q ■ tt^ 1 ■ q is a dimensionless scalar. 
We have found this form of the entropy for several classes 
of skewed trial solutions to Boltzmann's kinetic equation. 
Before we justify this ansatz and calculate the skewness 
contribution S(tp) to the entropy per particle in units of 
/cb for a particular class of trial functions, we give the 
derivatives of the entropy with respect to tt and q: 
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b(<p) tt 1 ■ qq ■ tt 1 



(3) 
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and 



SS 
5q 



— = —2nk-sb((p) 7r 1 ■ q, 



(4) 



with b((p) — —dS(ip)/d(p > 0. 

A class of 13 parameter functions. — We next derive the 
entropy @ for a natural class of approximate solutions 
to Boltzmann's kinetic equation. Our construction of 
approximate solutions to Boltzmann's kinetic equation 
is based on multivariate Gaussian distributions (see, for 
example, Sec. 2.1.3 of |13|). 



fx&(p) 



exp 



-(p-A).0- 1 -(p-A) 



(5) 



with a first-moment vector A, a positive-definite sym- 
metric second- moment tensor ©, and the normalization 
constant J\f = (2-7r) 3 / 2 (det 0) 1 / 2 . As a simpler alterna- 
tive to the Pearson- Type-IV distribution considered by 
Torrilhon we here propose to use the following class 
of skewed single-particle distribution functions, 

f{p) = n 1 - w ■ 1 ■ (p- A) 

x a((p-\)-&- 1 -(p-X))]f x& (p), (6) 

where a suitably chosen function a(x) regularizes the 
skewness introduced by w. A convenient choice of a(x) 
is discussed below. Note that / is nonnegative if ■& — 
w ■ © 1 • w < [pa(p 2 )]~ 2 for all p; it is hence important 
that there exists an upper bound for |pa(p 2 )| and, there- 
fore, that a(p 2 ) is nonpolynomial. In three dimensions, 
the scalar number density n, the vector A, the symmetric 
tensor 0, and the skewness vector w add up to a total 
of 13 parameters. These parameters are actually fields 
depending on the position r. The single-particle distri- 
bution functions © can be thought of as trial solutions 
to Boltzmann's kinetic equation or as an approximate 
invariant manifold 

A nice feature of our trial functions (JB]) is that all even 
moments coincide with those of the Gaussian distribution 
fx®- The odd moments can be reduced to more compli- 
cated Gaussian averages involving a(x). For the lowest 
three moments, we obtain the closed-form expressions 



M = n(X — ct4 w) 



© 



0(D) w, 0(tf) = (4a 4 - 2a 6 - 



a 4 ww, 



3ag — <i4 



(7) 



(8) 



with the numerical coefficients 
1 



2tt u-iy.U- 



p> a(p 2 )e-P 2 / 2 dp, (9) 



where fc!! is the product of all odd numbers from 1 to 
k. A convenient criterion for the choice of a{x) is that 
the numbers a, can be evaluated in closed form. For 




FIG. 1: Skewness contribution S{ip) to the entropy (continu- 
ous line) and linear approximation to it (dashed line). 



a{x) — exp(— ux/2), with a parameter u > 0, we have 
aj = (1 + u)~(i +1 ^ 2 and hence 04/06 = 1 + u. Once 
the relationship between the scalars $ = w ■ _1 • w and 



q ■ tt 



q is established, 



i) 



ip = 



1 



6(d) 2 , 



(10) 



the inversion of Eqs. ([T])-© to obtain the parameters of 
the trial solutions in terms of the moments is simple. 

By integrating —/In/ and using the relationship 
det7r = (1 — a 2 i!))det®, we confirm Eq. ([2]) and find 
the following skewness contribution to the entropy: 



1 



2tt jo 

2l 



zpa 



- ln[l — (zpa) 



1 \ 1 + zpa 
zpa J 1 — zpa 

I } p 2 e~P 2 / 2 dp -Un(l- a 2 (11) 



where z = yfd and a = a(p 2 ). For our example a(x) — 
exp(—ux/2) with the value u = 2/5, suggested to repro- 
duce the trial solutions of Grad's approach upon Taylor 
expansion of a(x), the corresponding entropy contribu- 
tion S(tp), as obtained by numerical integration and the 
relation (|10[) between ?? and tp, is represented by the con- 
tinuous line in Fig. [1] The entropy clearly favors the 
symmetric distributions characterized by (p — 0. The 
linear approximation S(ip) = —b(0) ip with 



4(0) = B5 



l + 2u 




(12) 



implying a quadratic dependence of S on q, is shown by 
the dashed line in Fig. [T] It provides a good approxi- 
mation over the entire range of allowed values of if. In 
the limit of small u, 6(0) is 1/20; for u — 2/5, we have 
6(0) ps 0.058. 

Evolution equations. — With an entropy at hand, we 
can now use our thermodynamic guiding principles to 
formulate evolution equations for our 13 moments. In 
addition to the usual continuity and momentum balance 
equations, 



dp 

m 



d_ 
dr 



(vp), 



dM 



d_ 
dr 



vM 



P 



(13) 
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where v = M/p is the velocity field, we have the fol- 
lowing evolution equation for the second moment tensor 

TV 

div d T 1 9 

— = -V ■ — TV - K ■ TV - TV • K — ■ [pQ) 

at or mp or 



1 



tr7r 



(14) 



where we have introduced the velocity gradient tensor 
with components Kjk = dvj/dr^. The reversible terms 
in the first line are exact, but we still need a closure 
approximation for the third moment tensor Q. From the 
general tensor structure, we expect 

Qijk = Q\{qiKjk + QjKik + qknj) - Q2 ii%q_k, (15) 

with suitable functions Q\ and Qi of tp (or, in view of 
Eq. (ITUl) . equivalently of d). For our trial functions ©, 
we obtain 



Qx = 



a 4 — ciq 



(Sag — 04) a\ 



(16) 



For •& = ip = 0, we find Qi = 1/5, and we suggest that 
2/5 is a realistic choice of Qi- Alternatively, Q\ and Qi 
could be obtained by matching simulation results or some 
established 13 moment equations. We here do not discuss 
further restrictions for the most general functional form 
of Qi and Qi that might result from the Hamiltonian 
character of reversible dynamics. The relaxation term in 
Eq. (frjj has the simplest possible form consistent with 
energy conservation in the collision process. 

We now turn to the evolution of q. From the perspec- 
tive of kinetic theory, q is a well-defined combination of 
second and third moments. From a phenomenological 
perspective, we only know that q is a covariant vector. 
A clear physical meaning is given to q through the closure 
(fT5"j) which relates q to the heat flux, 



P 

2rrv 



Q : 1. 



(17) 



We could, for example, fix the meaning of q by choosing 
Qi to be identically equal to 1/5. As an alternative, we 
could make q phenomenologically well-defined by choos- 
ing the form of the entropy contribution S(ip), for exam- 
ple, by choosing b to be constant (say, 1/20). As there 
is no straightforward phenomenological relationship be- 
tween third moments and entropy, we might even want 
to prescribe both Q\ an d b. As an evolution equation for 
the skewness vector q, we now propose 



dq 
~dt 
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-v ■ -rr-q 
or 
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k ■ q 
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q ■ tv 



d_ 
Or 



(pQ) 



1 



q tv 



Ampipb 
F 

m \ (p 



Tr ' 



(18) 



-11-22. w -i 



d 1 

-— tr7T D ■ q. 

dr 2t y 



All reversible terms in the first two lines of Eq. (p~8|) are 
dictated by our first thermodynamic guiding principle, 
that is, the conservation of the entropy ([2]) by reversible 
dynamics. The exact evolution equation for q involves 
divergences of the third and fourth moment tensors. The 
terms involving the divergence of pQ in Eq. (|18|) have 
exactly the same tensor structure as in the exact evolu- 
tion equation. The reversible term involving the scalar 
fitting parameter F represents the effect of fourth mo- 
ments. As this effect occurs only in the evolution of q 
but not of tv, it is highly restricted by our second guiding 
principle, the Hamiltonian character of reversible dynam- 
ics. The proposed term is the only one we could come 
up with in view of our two guiding principles. Whereas 
it does not match the full structure of the terms result- 
ing from fourth-moment closure, matching the hydrody- 
namic limit suggests F = 5/3. Finally, the relaxation 
term in Eq. (ITS)) involves a tensor D. 

From the evolution equations ([TBI . (IT4l . and (fT5|) we 
obtain the entropy balance equation 



ds 



d_ 
dr 



(vs) 



"67 



nk^b 



(tr7rtr7r _1 -9)H —q-rv^Dq, 

(19) 

where D = D — 1 + (tr7r/3)7r 1 and s is the entropy 
density implied by Eq. ([2]). For positive definite tv^ 1 ■ D, 
we have thus established an H theorem for the thermo- 
dynamically admissible 13 moment equations. From now 
on, we simply assume an isotropic tensor D = Dl. 

Hydrodynamic limit. — Within the Chapman-Enskog 
procedure for obtaining hydrodynamic equations in the 
limit of short relaxation times r, Eq. (Tl4"|) leads to 

tv - — 1 = -mk B TT{K + k t ), (20) 
3 

where T is the local equilibrium temperature implied by 
the trace of tv. Equation (|2U|) represents Newton's ex- 
pression for the stress tensor. For the discussion of heat 
flow, it is useful to look at the evolution equation for the 
scalar <p obtained by combining Eqs. (fT4|) and (TIB)) . 



dip dip 1 

— = —v ■ 

dt dr 2mpb 



--¥>• (2D 
r 



The Chapman-Enskog procedure then leads to 



ip 



d 



ik B TbD dr 



■3 



(22) 



After introducing this limiting form of ip into the last 
term of the entropy balance (|19l) , the proper combination 
of the thermal entropy flux and production arises. 

When Eq. (|22|) is used in the further analysis of 
Eq. (TIE)) , one can recognize that the term with prefac- 
tor F must become small. This observation implies that 
the temperature gradient becomes equal to its projection 
along q, that is, q must be proportional to the tempera- 
ture gradient. In other words, the direction of q is found 
by suppressing components perpendicular to the temper- 
ature gradient. 
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Summary and outlook. — The evolution equations (1131) , 
([T3]). and (fTS)) developed in this letter may be consid- 
ered as thermodynamically closed 13 moment equations 
for rarefied gas flow. They come with the nondecreasing 
entropy ([2]) and a physically appealing Hamiltonian for- 
mulation of reversible dynamics. Moreover, standard hy- 
drodynamics is recovered as a limiting case. As thermo- 
dynamics is designed exactly for that purpose, we expect 
well-behaved solutions in all physically meaningful appli- 
cations of these equations. However, a detailed study of 
solutions is beyond the scope of this fundamental devel- 
opment of thermodynamically admissible equations. 

From the perspective of thermodynamic modeling of 
skewness on the 13 moment level, the numbers or func- 
tions Qi, Q2, and F arising from the third and fourth 
moment closure and the entropy contribution S((p) [or, 
equivalently, its derivative b(ip) occurring in the evolution 
equation (|18[) ] may be considered as independent inputs. 
From the perspective of kinetic theory, these inputs (with 
a grain of salt in the case of F) are related to the function 
a(x) in the trial functions ([5]), namely by Eqs. © and 
(fTTI) . As long as there is no profoundly justified choice 
for a(x), however, one should feel free to fit all thermo- 
dynamic parameters or to use a simple set of parameters, 
for example Qi = 1/5, Q 2 = 2/5, F = 5/3, and b = 1/20. 



The dissipative parameters r and D can be obtained by 
matching our evolution equations with kinetic theory or 
experimental results. For Maxwell molecules, we expect 
L> = 4/3 

For most applications, the 13 moment equations need 
to be supplemented by boundary conditions. In the work 
[ij , the usefulness of thermodynamics in providing mean- 
ingful boundary conditions has been recognized. It may 
hence be worthwhile to point out that the GENERIC 
framework, employed here to obtain the bulk equations, 
has recently been supplemented by boundary thermody- 
namics [13, Gil- 

The treatment of rarefied gases on the level of 13 mo- 
ments, which is intermediate between kinetic theory and 
hydrodynamics, is a classical example of multiscale mod- 
eling. In addition to the thermodynamic exploration of 
this intermediate level, we have connected its parameters 
to kinetic theory. Our remarkably simple equations (|13|) , 
(|14"|). and (fTS)) show that thermodynamic principles, al- 
though strictly speaking never indispensable, offer enor- 
mously useful guidance. 
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